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Abstract 

Mixing patterns in large self-organizing networks, such as the Internet, the World Wide Web, social 
and biological networks are often characterized by degree-degree dependencies between neighbouring 
nodes. In this paper we propose a new way of measuring degree-degree dependencies. One of 
the problems with the commonly used assortativity coefficient is that in disassortative networks 
its magnitude decreases with the network size. We mathematically explain this phenomenon and 
validate the results on synthetic graphs and real-world network data. As an alternative, we suggest 
to use rank correlation measures such as Spearman's rho. Our experiments convincingly show that 
Spearman's rho produces consistent values in graphs of different sizes but similar structure, and 
it is able to reveal strong (positive or negative) dependencies in large graphs. In particular, we 
discover much stronger negative degree-degree dependencies in Web graphs than was previously 
thought. Rank correlations allow us to compare the assortativity of networks of different sizes, which 
is impossible with the assortativity coefficient due to its genuine dependence on the network size. We 
conclude that rank correlations provide a suitable and informative method for uncovering network 
mixing patterns. 
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I. INTRODUCTION 



This paper proposes a new way of measuring mixing patterns in large self-organizing 
networks, such as the Internet, the World Wide Web, social and biological networks. Most 
of these real- world networks are scale-free, i.e., their degree distribution has huge variability 
and closely follows a power law (the fraction of nodes with degree k is roughly proportional 
to k~^~^, 7 > 0). We study correlations between degrees of two nodes connected by an edge. 
This problem, first posed in [H, Hj, has received vast attention in the networks literature, in 
particular in physics, sociology, biology and computer science. We show however, analytically 
and on the data, that the presence of power laws makes currently used measures inadequate 
for comparison of mixing patterns in networks of different sizes, and provide an alternative 
that is free from this disadvantage. 

Adequate measuring and comparison of degree-degree correlations is important because 
mixing patterns define many of the network's properties. For instance, the Internet topology 
is not sufficiently specified by the degree distribution; the negative degree-degree correlations 
in the Internet graph have a great infiuence on the robustness to failures j3|, efficiency of 
Internet protocols [j], as well as distances and betweenness j5|. This is totally different from 
the mixing patterns in networks of bank transactions [6| where the core of 25 most important 
banks is entirely connected. The correlation between in- and out-degree of tasks plays and 
important role in the dynamics of production and development systems 0]. Mixing patterns 
affect epidemic spread p, 0] and Web ranking [10] . 

In his seminal papers, Newman [ll, |2| proposed to measure degree-degree correlations using 
the assortativity coefficient, which is, in fact, an empirical estimate of the Pearson's correlation 
coefficient between the degrees at either ends of a random edge. A network is assortative 
when neighbouring nodes are likely to have a similar number of connections. In dis assortative 
networks, high-degree nodes mostly have neighbours with small number of connections. The 
empirical data in [ll. Table I] suggest that social networks tend to be assortative (which is 
indicated by the positive assortativity coefficient), while technological and biological networks 
tend to be disassortative. 

In [H, Table I], it is striking that larger disassortative networks typically have an assorta- 
tivity coefficient that is closer to and therefore appear to have approximately uncorrelated 
degrees across edges. Similar conclusions can be drawn from P, Table II] . In recent literature 



111 . |12| the issue was raised that the Pearson's correlation coefficient in scale-free networks 
decreases with the network size. In this paper we demonstrate analytically and on the data 
that in all scale- free disassortative networks with a realistic value of the power-law exponent, 
the assortativity coefficient decreases in magnitude with the size of the graph. In assortative 
networks, on the other hand, the assortativity coefficient can show two types of behaviour. It 
either decreases with graph size, or it shows a considerable dispersion in values, even if large 
networks are constructed by the same mechanism. 

We suggest an alternative solution based on the classical Spearman's rho measure [l3| that 
is the correlation coefficient computed on the ranks of degrees. The huge advantage of such 
dependency measures is that they work well independently of the degree distribution, while 
the assortativity coefficient, despite the fact that it is always in [—1, 1], suffers from a strong 
dependence on the extreme values of the degrees. The usefuUness of the rank correlation 
approach to discover dependencies in skewed distributions has already been postulated in the 
1936 paper by H. Hotelling and M.R. Pabst [13] : 'Certainly where there is complete absence 
of knowledge of the form of the bivariate distribution, and especially if it is believed not to be 
normal, the rank correlation coefficient is to be strongly recommended as a means of testing 
the existence of relationship. ' 

We compute Spearman's rho on artificially generated random graphs and on real data 
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from web and social networks. Our results agree with concerning the presence of positive 
or negative correlations, but Spearman's rho has two important advantages: (1) it is able to 
reveal strong disassortativity in large networks; (2) it produces consistent values on the graphs 
created by the same mechanism, e.g. on preferential attachment graphs of different sizes. 
Thus, Spearman's rho correctly and consistently captures the underlying connection patterns 
and tendencies. We conclude that when networks are large, or two networks of difference sizes 
must be compared (e.g. in web crawls or social networks from different countries). Spearman's 
rho is a preferred method for measuring and comparing degree-degree correlations. 

The closing section discusses further challenges in the evaluation of network mixing pat- 
terns. 



II. NO DISASSORTATIVE SCALE-FREE RANDOM GRAPH SEQUENCES 

In this section we present a simple analytical argument that in disassortative networks 
the assortativity coefficient always decreases in magnitude with the size of the graph. Formal 
proofs can be found in 

Assortativity in networks is usually measured using the assortativity coefficient, which is 
in fact a statistical estimator of a Pearson's correlation coefficient for the degrees on the two 
ends of an arbitrary edge in a graph. Let G = {V, E) be a graph with vertex set V, where 
|\^| = n denotes the size of the network, and edge set E. The assortativity coefficient of G is 
equal to (see, e.g., (4)]) 
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where the sum is over directed edges of G, i.e., ij and ji are two distinct edges, and is the 
degree of vertex i. We compute that 



^ 5: |(rf. + d,) = E wii: -M + d]) = ^ E ^^ 

Thus, pn can be written as 

Sije-B didj ~ |E| ( Siev df^ 
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In practice, all quantities in (1II.2P are finite, and p„ can always be computed. However, 
since many real-life networks are very large, a relevant question is how p„ behaves when n 
becomes large. 

In the literature, many exarnples are reported of real-world networks where the degree 
distribution obeys a power law 17, 3]. In particular, for scale-free networks, the observed 
proportion of vertices of degree k is close to /(fc ) = Cq/c"'^"^, and most values of 7 found 
in real-world networks are in (1,3), see e.g., [171, Table I] or 18|, Table I]. For p < 7, let 
/ip = ^fcfc^/(fc), and note that the series diverges if p > 7; let a ~ 6 denote that ajh — )■ 1. 
Then we can expect that, as n grows large, 

\E\ = ^^di^ pxfi, E'^?~/^p^' P < 7> 
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while maxjev' di is of the order n^^'^. As a direct consequence, 

cn < \E\ < Cn, (II.3) 

cn^/^ < maxrf, < Cn^/^, (II.4) 

i€[n] 

^^max{p/7,i} < ^ rfP < c-n-^^^^P/^'i}, p = 2, 3, (II.5) 

ie[n] 

for 7 G (1,3) and some constants < c < C < oo. We emphasize that conditions f lII.3P - 
( ]II.5|) are very general and hold for any scale-free network of growing size, independently of 



its mixing patterns. From ( 1II.2P we simply write 
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and notice that 



whereas 



> (maxrfi)^ > c^n^/^. 



Since 7 G (1, 3) we have max{4/7 — 1, 1} < 3/7, so that 

Xliey df 
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—7- 00 as t —7- 00. 



Hence, the lower bound p~ is of the order 72,™^^{i/t~1'1^3/7}_ jg easy to check that if 
7 G (1,3), then p~ converges to zero when the graph size increases. This means that any limit 
point of the assortativity coefficients p„ is non-negative. Note also that p~ is defined by the 
degree sequence, and it does not depend on the mixing pattern at all. We conclude that by 
looking only at the value of pn one cannot discover even very strong disassortativity in large 
scale-free graphs. We will confirm this finding in Section |IV] on artificially generated random 
graphs, and in Section |V] on real-world networks. 

We note that if 7 > 3, then all terms in f lll.ip converge to a number, and p„ does not 
scale with the network size. In practice this means that the dependence of p„ on the graph 
size is observed when node degrees have a broad distribution, and this range increases when 
the network gets bigger. This is the case in most real-life networks and models for them, as 
is e.g. obviously the case for preferential attachment models. 

We further notice that f lII.3p -( 1lL5|) imply that 

^ didj < (^maxc/i) J2 = maxd,(^J2 dtj < C2ni/^+'^^"^2/7,i}_ (ng) 

Mathematically, an interesting case is when '^^j^^ d^dj and "^^^y df are of the same order of 
magnitude. Then the network is assortative but, formally, p„ converges to a random variable. 
In practice this means that p„ can result in very different values on two very large graphs 
constructed by the same mechanism. We will give such an example in Section IIVI 
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III. RANK CORRELATIONS 



We propose an alternative measure for the degree- degree dependencies, based on the 
rank correlations. For two-dimensional data ((Xj, Fj))"^^, let rf and be the rank of an 
observation Xi and Yi, respectively, when the sample values (Xj)"^^ and are arranged in 

a descending order. The rank correlation measures evaluate statistical dependences on the data 
{{rf, rY))f=i, rather than on the original data {{Xi, Yi))^^^. Rank transformation is convenient, 
in particular because {rf ) and {rj) are samples from the same uniform distribution, which 
implies many nice mathematical properties. 



The statistical correlation coefficient for the rank is known as Spearman's rho [13 



.rank _ ELll^'' " + ^) ^1 " + l)/2) 



rank ^ A^z=l\ i \ ' Jl J\ t \ ' Jl J /jjj 

^y::=m - (^+ i)/2)^E:(rf - {n+i)/2f 

The mathematical properties of the Spearman's rho have been extensively investigated. In 
particular, if {{Xi, l^i))r=i consists of independent realizations of (X, Y), and the joint distribu- 
tion function of X and Y is different iable, then p™'^ is a consistent statistical estimator, and 
its standard deviation is of the order 1 / y/n independently of the exact form of the underlying 



distributions, see e.g. [19 



For a graph G of size n, we propose to compute p^^^^ using (lIII.ip as follows. We define 
the random variables X and Y as the degrees on two ends of a random undirected edge in a 
graph (that is, when rank correlations are computed, ij and ji represent the same edge). For 
each edge, when the observed degrees are a and 6, we assign [X = a,Y = h] oi [X = h,Y = a] 
with probability 1/2. Many values of X and Y will be the same making their rank ambiguous. 
We resolve this by we adding independent uniformly distributed random variables on [0, 1] to 
each value of X and Y . In the setting when the realisations {Xi, Yi) are independent, this 
way of resolving ties preserves the original value of the Spearman's rho on the population. 



see e.g. 01. We refer to [2l| for a general treatment of rank correlations for non-continuous 
distributions. 

In the remainder of the paper we will demonstrate that the measure pj^f"'^ gives consistent 
results for different n, and it is able to reveal strong negative degree-degree correlations in 
large networks. 



IV. RANDOM GRAPH DATA 



We consider four random graph models to highlight our results. 

The configuration model. The configuration model was invente d by Bollobas in 23], in- 



spired by ^23|]. It was popularized by Newman, Strogatz and Watts [24|, who realized that it 
is a useful and simple model for real-world networks. In the configurations model a node i 
has a given number di of half-edges, with Cn = ^jgy di assumed to be even. Each half-edge is 
connected to a randomly chosen other half-edge to form an edge in the graph. We chose 7 = 2, 
thus, the maximum degree is of the order n^^^, which corresponds to the case of uncorrelated 
random networks, such that the probability that two vertices are directly connected is close 
to didj/in 25|, |26|. Although self-loops and multiple edges can occur these become rare as 
n oo, see e.g. [27| or 28|. In simulations, we collapse multiple edges to a single edge, and 
remove self-loops. This changes the degree distribution slightly, and intuitively should yield 
negative dependencies. In Figure [U^a) we observe that, on average, p„ and pl^^^ are indeed 
negative in smaller networks but then they converge to zero showing that the degrees on two 
ends of a random edge are uncorrelated. 
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Configuration model with intermediate vertices. In order to construct a strongly disassor- 
tative graph, we first generate a configuration model as described above, and then we replace 
every edge by two edges that meet at a middle vertex. In this model, there are n + inf^ 
vertices and 2£„ edges (recall that ij and ji are two different edges). Now, if E, V, and di, 
i = 1, . . . ,n denote, respectively, the edge set, the vertex set, and the degrees of the original 
configuration model, then in the model with intermediate edges the assortativity coefficient is 
as follows: 



Pn 



2 • 



When 7 < 3 we have ps = oo, and thus p„ — > as n — )■ cxd. Furthermore, the lower bound p~ 
also converges to zero as n grows. It is clear that this particular random graph, of any size, 
is equally and strongly disassortative, however, p„ fails to capture this. In Figure [T]^b) it is 
clearly seen that both p„ and p~ quickly decrease in magnitude as n grows. It is striking that 
^rank gj^Qws a totally different and very appropriate behavior. Its values remain around —0.75 
identifying the strong negative dependencies, and the dispersion across different realizations 
of the graph decreases as n — t- cxo. 

Preferential attachment model. We consider the basic version of the undirected prefer- 
ential attachment model (PAM), where each new vertex adds only one edge to the network. 



connecting to the existing nodes with probability proportional to their degrees [15|. In this 
case, it is well known that 7 = 2 (see e.g. ^]). Newman [1| noticed the counterintuitive fact 
that the Preferential Attachment graph has asymptotically neutral mixing, p„ — > as n — )■ 00. 



This phenomenon has been studied in detail by Dorogovtsev et al. and it can be clearly 
observed in Figure [H^c). The reason for this behavior is not the genuine neutral mixing in 
the PAM but rather the unnatural dependence of p„ on the graph size. Indeed, we see that 
PAMs of small sizes have p„ < 0, and then the magnitude of p„ decreases with the graph size. 
Again, Spearman's rho consistently shows that the degrees are negatively dependent. This 
can be understood by noting that the majority of edges of vertices with high degrees, which 
are old vertices, come from vertices which are added late in the graph growth process and 
thus have small degree. On the other hand, by the growth mechanism of the PAM, vertices 
with low degree are more likely to be connected to vertices having high degree, which indeed 
suggests negative degree-degree dependencies. 

A collection of complete bipartite graphs. We next present an example where the assor- 
tativity coefficient has a nonvanishing dispersion. Take {{Xi,Yi))^^i to be a sample of inde- 
pendent realizations of the vector (X, Y). We assume that X = bUi + bU2 and Y = bUi + af/2, 
where 6 > 0, a > 1, and Ui,U2 are independent identically distributed (i.i.d.) random vari- 
ables with power law tail, and tail exponent 7. Then, for z = 1, . . . ,n, we create a complete 
bipartite graph of Xi and Yi vertices, respectively. These n complete bipartite graphs are not 
connected to one another. We denote such a collection of n bipartite graphs by Gn- This 
is an extreme scenario of a network consisting of highly connected clusters of different size. 
Such networks can serve as models for physical human contacts and are used in epidemic 
modelling [9|. 



The graph Gn has \V\ = ^"=i(-^i + Yj) vertices and \E\ = 2 Y^^=i ^lYi edges. Further, 

n n 

ieV i=l ij&E 1=1 

Assume that P(f/j > x) = cqx'"' , where cq > 0, x > xq, and 7 G (3, 4), so that E[?7^] < 00, but 
E[[/^] = 00. As a result, \E\/n 2E[Xr] < cx) and i Y.ievdf ^ E[Xr(X + Y)] < 00. 
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Further, 



N 

i=l 



where Zi and Z2 and two independent stable distributions with parameter 7/4. As a result, 



d 2a'Zi + 2^2 



Pn — ^ 7 — , , , as n ^ 00, 

[a + a'^jZi + 2Z2 



which is a proper random variable taking values in (2a/ (1 + a?), 1), see 16| for detailed proof. 

Note that in this model there is a genuine dependence between the correlation measure 
and the graph size. Indeed, if n = 1 then the assortativity coefficient equals —1 because 
nodes with larger degrees are connected to nodes with smaller degrees. However, when the 
graph size grows, the positive linear dependence between X and Y starts dominating, thus, 
larger graphs of this structure are strongly assortative. While the example we present is quite 
special, we believe that the effect described is rather general. 

In Figure [It^d) we again see that p^^^^ captures the relation faster and gives consistent 
results with decreasing dispersion. On a contrary, p„ has a persistent dispersion in its values, 
and we know from the result above that this dispersion will not vanish as — )■ 00. In the 
limit, pn has a non-zero density on (0.8, 1). However, the convergence is too slow to observe 
it at n = 100, 000, because the vanishing terms are of the order n~^f^ ^ which is only 77,^1/31 
our example. 



V. WEB SAMPLES AND SOCIAL NETWORKS 

We computed p„, p^^^^ and p^ on several Web samples (disassortative networks) and 
social network samples (assortative networks). We used the compressed graph data from the 
Laboratory of Web Algorithms (LAW) at the Universita degli studi di Milano js^, 3l[. We 



used the bvgraph MATLAB package [32]. The stanford-cs database [33|] is a 2001 crawl that 
includes all pages in the cs.stanford.edu domain. In datasets (iv), (vii), (viii) we evaluate pn, 
p^^^^ and p~ over 1000 random edges, and present the average over 10 such evaluations (in 10 
samples of 1000 edges, the observed dispersion of the results was small). 

The results are presented in Table [B We clearly see that the assortativity coefficient p„ and 
Spearman's p^^'^'^ always agree about whether dependencies are positive or negative. They also 
agree in magnitude of correlations when graph size is small or the lower bound p~ is sufficiently 
far from zero. However, p„ is not consistent for graphs of similar structure but different sizes. 
This is especially apparent on the two .uk crawls (iii) and (iv). Here p„ is significantly smaller 
in magnitude on a larger crawl. Intuitively, mixing patterns should not depend on the crawl 
size. This is indeed confirmed by the value of Spearman's rho, which consistently shows strong 
negative correlations in both crawls. We could not observe a similar phenomenon so sharply 
in (vi) and (vii), probably because a larger co-authorship network incorporates articles from 
different areas of science, and the culture of scientific collaborations can vary greatly from one 
research field to another. 

We also notice that, as predicted by our results, the assortativity coefficient tends to take 
smaller values than p^^^'' if p~ is small in magnitude. This is clearly seen in the data sets 
(ii), (iv) and (v). Again, (ii) and (iv) are the largest among the analyzed web crawls. 

The observed behaviour of the assortativity coefficient is explained by the above stated 
results that p„ is infiuenced greatly by the large dispersion in the degree values. The latter 
increases with graph size because of the scale-free phenomenon. As a result, p„ becomes 
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FIG. 1. (Color online) Scatter plots for samples of 20 graphs. For each size we plot the 20 realizations 
of pn (blue asterisks) and p^n'^^ (^^"^ diamonds) in random graphs of different sizes. Solid lines connect 
the averages of the samples. In (c), (d) the circles connected by the solid line are the averages of 
p~ in the samples, (a) Configuration model, F{d > x) = x~'^, x > 1. (b) Configuration model with 
intermediate vertices, (c) Preferential attachment model, (d) A collection of bi-partite graphs, where 
h = 1/2, a = 2, and U has a generalized Pareto distribution P([/ > x) = ((2.1 + x)/3.1)~^'"'^, x > 1. 



nr 


Dataset 


Description 


# nodes 


# edges 


max degree 


Pn 


„iank 
rn 


Pn 


(i) 


stanford-cs 


web domain 


9,914 


54,854 


340 


-.1656 


-.1627 


-.4648 


(ii) 


cu-2005 


.eu web domain 


862,664 


5,477,938 


68,963 


-.0562 


-.2525 


-.0670 


(iii) 


uk@100,000 


.uk web crawl 


100,000 


5,559,150 


55,252 


-.6536 


-.5676 


-1.117 


(iv) 


uk@l,000,000 


.uk web crawl 


1,000,000 


77,123,940 


403,441 


-.0831 


-.5620 


-.0854 


(v) 


enron 


e-mail exchange 


69,244 


506,898 


1,634 


-.1599 


-.6827 


-.1932 


(vi) 


dblp-2010 


co-authorship 


326,186 


1,615,400 


238 


.3018 


.2604 


-.7736 


(vii) 


dblp-2011 


co-authorship 


986,324 


6,707,236 


979 


.0842 


.1351 


-.2963 


(viii) 


hollywood-2009 


co-starring 


1,139,905 


113,891,327 


11,468 


.3446 


.4689 


-0.6737 



TABLE I. (i)-(iv) Web crawls: nodes are web pages, and an (undirected) edge means that there 
is a hyperlink from one of the two pages to another; (iii),(iv) are breadth-first crawls around one 
page, (v) e-mail exchange by Enron employees (mostly part of the senior management): node are 
employees, and an edge means that an e-mail message was sent from one of the two employees to 
another, (vi), (vii) scientific collaboration networks extracted from the DBLP bibliography service: 
each vertex represents a scientist and an edge means a co-authorship of at least one article, (viii) 
vertices are actors, and two actors are connected by an edge if they appeared in the same movie. 
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smaller in magnitude, which makes it impossible to compare graphs of different sizes. In 
contrast, the ranks of the degrees are drawn from a uniform distribution on [0, 1], scaled by 
the factor n. Clearly, when a correlation coefficient is computed, the scaling factor cancels, 
and therefore Spearman's rho provides consistent results in the graphs of different sizes. 



VI. DISCUSSION 

The assortativity coefficient pn proposed in [l[ |2[ has been the first dependency measure 
introduced to describe degree-degree correlations in networks. The assortativity coefficient 
has provided many interesting insights. It has been successfully used for comparison of de- 



pendencies in graphs with the same degree sequences [3J, |35[, and to generate graphs with 



given degrees and desired mixing patterns [36|. An important drawback of pn is its depen- 



dence on the network size n. It has been noticed by many authors, and shown in this paper 
for disassortative networks, that p„ converges to zero as n grows. In particular, the decay with 
network size of the assortativity coefficient p„ implies that it cannot be used for comparing 
dependencies in networks of different sizes. Therefore, it prohibits the investigation whether 
growing networks become more or less assort it at ive over time. 

This paper suggests to use rank-correlation measures such as Spearman's rho. Our ex- 
periments convincingly show that Spearman's rho does not suffer from the size-dependence 
deficiency. In networks of different sizes but similar structure. Spearman's rho yields consistent 
results, and it is able to reveal strong (positive or negative) correlations in large networks. We 
conclude that rank correlations are a suitable and informative method for uncovering network 
mixing patterns. 

For the correct interpretation of degree-degree dependencies, it is important to realise 
that positive or negative correlations can be pre-defined by the degree sequence itself. For 
instance, there is only one simple graphs with degrees (3,1,1,1), and the result p^ = — 1 
is not informative in this case. It has been discussed in the literature that, conditioned 
on not having self-loops and multiple edges, random networks with given degrees exhibit 
disassortative patterns j25|, [ss], S^l , also called structural correlations. In order to filter out 



the structural correlations, one needs to compare the real- world networks to their null-models - 
graphs with the same degree sequences but random connections. This null-model is a uniform 
simple random graph with the same degree sequence. Here a network is called simple when 
it has no self-loops nor multiple edges. Such a graph can be obtained by randomly pairing 
half-edges, as in Section IIV[ and taking the first realization that is simple. This is especially 
problematic when (maxj c/j)^ > \E\, which is the case in many examples, since then one needs 



a prohibitingly large number of attempts before a simple graph is generated [28|, |38 



A widely accepted method for constructing a null-model, is the random rewiring of the 



connections in a given graph |34j . |35| . The disadvantage is the unknown running time before 



a graph is produced that is close enough to being uniform. Recent work [39[ presents a 
sequential algorithm, where, at each step, the remaining unconnected edges maintain the 
ability to generate a simple graph. This method always produces the desired outcome but its 
worst-case running time 0{n'^^^di) is infeasible for large networks. The recently introduced 
grand-canonical model [ioj computes the probability of connection between two nodes in a 
maximum entropy graph with given degree sequence, and enables the evaluation of many 
characteristics of the graph. To the best of our knowledge, efficient implementation of this 
method for large networks has not been developed yet. 

Constructing a null-model and filtering out the structural correlations in large networks is 
an interesting and demanding computational task that is beyond the scope of this paper. We 
believe that structural correlations will affect p„ to a larger extent than the rank correlation 
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prank because it is usually the nodes with largest degrees that produce self-loops and multiple 
edges, and thus the relative contribution of these edges in the cross-products will be larger for 
Pn than for pj^^"*^. This conjecture requires a further investigation. 

We conclude by stating that rank correlation measures deserve to become a standard 
tool in the analysis of complex networks. The use of rank correlation measures has become 
common ground in the area of statistics for analysing heavy-tailed data. We hope to have 
provided a sufficient evidence that this method is preferred for analysing network data with 
heavy-tailed degrees as well. 
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